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Introduction 
Introduction to the Direction of Arrival Estimation beamforming ELEC 301 
project 


Introduction 


Beamforming is a type of signal processing which uses sensor arrays to 
either transmit or receive signals in a certain direction. This directionality is 
accomplished by manipulating the phases at each sensor in the array so that 
constructive interference occurs at favored angles and destructive 
interference occurs elsewhere. Beamforming is most often used with sound 
or radio waves and it has numerous applications in areas such as RADAR, 
SONAR, wireless communications, and acoustics. 


Direction of Arrival Estimation 


Our specific application of beamforming is direction of arrival estimation. 
This process is described in detail in the module of that name, but it can be 
roughly summed up in a sentence as using a sensor array to determine at 
what angle, in relation to the center of the array, a propagating source exists. 
For our simulation we used a linear array of sensors along the x-axis and 
defined the phase center of the array as the point (0, 0). If there were an odd 
number of sensors the center sensor was on the origin, and if there were an 
even number the center was between the two central sensors. Here is a way 
to visualize the setup of our MATLAB simulation: 


Signal source 


ceca 


Sensor Array 


The goal is to find the angle @. It may be helpful to keep this image in mind 
as the techniques behind beamforming and direction of arrival are 
explained. 


Beamforming 
A discussion of the conventional frequency-domain beamforming algorithm 


Beamforming 


The conventional frequency domain beamforming algorithm involves three 
components: a steering vectore, a weighting matrix W, and a vector Y 
which contains the frequency information about the signal received at each 
sensor. The steering vector is given its own section, and I will go into brief 
detail about the other two here. 


Sensor Data 


The column vector Y is actually a collection of short-time Fourier 
Transforms. The signal received at each sensor at a given time t is 
multiplied by a window of duration D and then the Fourier Transform of 
that product is taken. Mathematically, this can be expressed as 


Y,(t,«) = | S(t — Dy, (de*tdr 


This transform effectively represents the frequency components of the 
signal during various blocks of time, and is useful if the signal itself is finite 
and/or not uniform over all time. When performing a short-time Fourier 
Transform, it is best to pick a window that has a Fourier Transform which 
looks like a single narrow spike so that its width is narrow compared to the 
signal’s spectral variations — this allows the short-time Fourier Transform of 
a signal at one instant in time to best approximate the signal’s spectrum. 
The vector Y should actually be written as Y(t, @) because it is dependent 
both on time and the frequency content of the signal, but often that notation 
is omitted for neatness. Each element Yn is the collection of short-time 
Fourier Transforms at sensor n, thus giving Y the form 


¥,(t @) 
Y(t,@) = ove 
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The Weighting Matrix 


The matrix W is diagonal; the only elements which are not zero are the 
sensor weights along its main diagonal. The actual content of those weights 
can vary by application. In a simple linear array like the one used in our 
particular example, all of the weights are simply one. However, in more 
interesting array configurations such as an array with sensors that were 
placed randomly or real-world sensors that have inherent noise components, 
the weighting matrix can become much more important. 


Beamformer Output 


The Fourier Transform of the beamformer output can be written in the 
following way: 
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This is the inner product of e’ and WY, which is nothing more than a 
column vector of weighted short-time Fourier transforms. In order to get the 
beamformer output in the time domain, all that is necessary is taking an 
inverse short-time Fourier Transform of e’WY‘(t, @). 


Here is a visual representation of the entire beamforming process: 
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Steering Vector 
A brief overview of what a steering vector is and why it is important to 
beamforming 


The Steering Vector 


The steering vector e represents the set of phase delays experienced by a 
plane wave as it reaches each element in an array of sensors. The 
mathematical representation of it is 


exp{—jk® * Tal 


e= wens 

exp{—j 3 Xy-1} 
Where j is the imaginary unit, k, is the wave number at the nth sensor, and 
X, is the position vector of the nth sensor in (x, y, z) coordinates. If the 
sensor array is two-dimensional, all z coordinates can be treated as 0 and 
omitted from the position vector in calculations. The wave number of a 
plane wave is its spatial frequency, or how many waves exist over a 
specified distance. It is a function of the frequency of the propagating wave, 
the material through which it is travelling, and its original location or 
direction of propagation. For our specific application, we actually use the 
angular wave number, which is written as (27/A) and indicates the number 
of wavelengths per 27 unit distance. 


The role of the steering vector e in beamforming is to define the assumed 
propagation direction of the source’s signal. In other words, it is the part of 
the beamforming process that chooses the direction in which the sensor 
array creates an output beam. Because e is designed to cancel the phase 
shifts at each sensor in the array and the wave number k must, by definition, 
assume a propagation direction, the steering vector has the effect of moving 
the each sensor’s output to a particular direction, thus focusing the beam. 


Matched Filter 
A description of matched filters, which are a component of our direction of 
arrival estimation algorithm. 


Matched Filtering 


Matched filters underlie many signal classification systems in the real 
world, allowing systems to determine the likelihood of a signal being 
present in a noisy environment. They take the form a cross-correlation or 
convolution between the received signal and a desired signal, the presence 
(or absence) of which is to be ascertained. This cross-correlation may then 
be compared to a threshold, comprising a decision rule to test for the 
presence of the known signal. 


Time-Domain Implementation 


Within the time domain, matched filters take the form of a convolution of 
the received (x[n]) and desired (s[n]) signals. 
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The output of this convolution will have a distinct peak when the reference 
and received signals are most closely correlated, and will have smaller 
fluctuations when they are uncorrelated or weakly correlated. Therefore, the 
time domain output of this convolution may simply be compared to a 
threshold. For samples where y[n] exceeds the threshold, we assume that 
x[n] contains our reference signal, and samples where y[n] falls below the 
threshold, x[n] does not contain our reference. 


Frequency-Domain Implementation 


Through the properties of the Fourier Transform, we can transform the 
above convolution into frequency-domain multiplication, giving rise to the 
frequency domain implementation of the matched filter. 


Y[k] = S[k] « X[k] 


While this implementation is obviously more computationally efficient if a 
frequency-domain copy of s[n] and x[n] already exists (as in the frequency- 
domain beamforming implementation), there is an obvious disadvantage for 
naive implementations. If one simply took the DFT over all incoming 
samples to yield X[k], then no localization of the reference signal in time 
would be possible. That is, a decision could possibly be made about the 
existence or absence of the reference within the overall signal, but it could 
not be localized to a specific time within that interval. 


The solution to this falls out of the Short Time Fourier Transform (STFT). 
If, instead of computing the DFT over all samples of the input, it is 
computed over a windowed subset, then we can still apply the above 
matched filter method while maintaining the ability to localize the signal’s 
occurrence in time. The new representation of this is as follows for non- 
time-varying s[n]: 


Y [k,n] = S[k] « X[k,n] 


The output of the frequency-domain matched filter can then be either 
transformed back into the time domain with the inverse STFT to be used 
with the existing per-sample decision rule, or can be used as the input to a 
frequency domain interpretation of the decision rule. The latter is obviously 
more computationally efficient, and is, at its simplest, the L? norm of Y[k, 
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P[n] may then be compared to a threshold to determine the presence or 
absence of the reference signal within the STFT slice n. 


Direction of Arrival Estimation 
A description of what direction of arrival is and a discussion of our specific 
algorithm for determining it 


Direction of Arrival Estimation 


Given the basic foundation in beamforming the other modules in this 
section have established, we would like to extend the usefulness of the 
technique into new domains. One application of interest outside of simply 
creating a directional antenna is to identify the direction from which a 
known source is broadcasting, called Direction of Arrival Estimation. 


Basic Premise: 


In all cases, direction of arrival estimation relies upon an a priori estimate 
of the signal’s arrival direction. If no estimate can be made, the array’s 
response is swept around the unit circle, and a decision rule is applied to its 
output. Locations where this decision rule is satisfied are assumed to 
correspond to locations of the known source. If an estimate exists, more 
sophisticated adaptive processing techniques may be applied to determine 
the source’s location with greater precision. 


Specific Implementation: 


To simplify the process of direction estimation, we make a few assumptions 
about the configuration of our sensor network and its environment. The first 
assumption we make is that our source is sufficiently distant from our 
sensor network for the received wavefront to be considered linear. This first 
assumption is called the Far Field Approximation, and simplifies the 
computation of steering and weighting vectors. The second assumption we 
make pertains to the noise environment the sensor operates within. We 
assume that noise (if present) is ideal Gaussian noise, simplifying the 
structure of our signal detection segment. 


These assumptions having been made, the detection algorithm is quite 
simple. We step through angles over which the array can look, feeding an 
input from a known location. The output of the array is then fed through a 
matched filter, which has the effect of computing the cross-correlation 
between the signal received by the array and the known signal of the 
transmitter. The peak value of that correlation is then stored with its 
corresponding angle. Upon completing the sweep over angles, the angle 
with the highest peak correlation is taken to be the most likely location of 
the signal source. This algorithm is described below in pseudocode: 


Pseudocode: 


Function AngleEstimation{ 
Angle <- 0 
Angles <- [] //empty set of angles 
Magnitudes <- [] //empty set of peak 
amplitudes 
While Angle < 2*pi { 
Output <- Beamformer(Angle, input) 
Output <- Matched_Filter (Output) 
Angles <- concatenate(Angles, Angle) 
Magnitudes <- concat(Magnitudes, 
max (Output) ) 
Angle <- Angle + pi/32 //increment 
to next angle 
} 
SourceAngle <- 
Angles[index(max(Magnitudes ) ) ] 
Return SourceAngle; 
; 


Conclusion 
A brief conclusion of our project and an overview of more advanced 
beamforming techniques 


Conclusion 


Beamforming is a sizeable and complex sub-study of digital signal 
processing. While this module has covered the fundamentals of a basic 
frequency domain implementation, more sophisticated techniques certainly 
exist. As such, this segment will serve less as a synopsis of the previous 
sections and more an overview of more advanced techniques in the field. 


Synopsis 


Beamforming is, at its core, the art of combining multiple omnidirectional 
sensors into a directional whole with superior characteristics to purpose- 
built directional sensors. This feat is performed by delaying the sensors’ 
outputs such that they constructively interfere for signals arriving from a 
specific direction, and may be done in either the time or frequency domain. 
While time domain implementations may be simpler for some array 
topologies, frequency-domain implementations are, by and large, much 
more efficient due to their ability to create fractional delays (a 
computationally intensive upsample-delay-downsample process in time 
domain). 


These systems underlie many technologies that we take for granted, from 

synthetic aperture RADAR used for such disparate applications as weather 
radar and ICBM tracking, to next-generation cellular technologies seeking 
to improve connectivity by steering cell signals towards active cell towers. 


Advanced Beamforming 


While the basic techniques described previously yield acceptable results 
(particularly for our purposes), adaptive processing can be used to far 


increase the performance of an array. These techniques largely involve a 
weighting vector w which was not used in our conventional beamforming 
algorithm. As this weighting vector may be thought of as a windowing 
function applied to the spatial characteristics of the array, it is clear that 
careful choice of a weighting function could improve parameters of the 
beamformer, such as its mainlobe width. 


The first method to perform this width reduction is called Minimum 
Variance beamforming, which can also be known as Capon’s method. This 
algorithm, as its name implies, seeks to find a weighting vector w such that 
the expectation value of |w’y|* is minimized. This is clearly the equivalent 
of minimizing the variance of the signal in zero-mean systems. Computing 
the value of w depends on the propagation direction and the spatiospectral 
covariance of all the sensors, but the technical details of that computation 
are outside the scope of this summary. Upon applying this weighting vector, 
the mainlobe width of the array’s output will be much smaller than that 
generated by the conventional beamforming algorithm. 


Even more advanced techniques, such as Linear Prediction, can yield even 
better results for the beamformer’s mainlobe width. At a high level, Linear 
Predication involves expressing the output of one particular sensor as a 
weighted linear combination of all the others. Again, the mathematical 
details of this method are beyond the scope of this module. 


Further adaptive methods allow for the aggregation of multiple arrays to 
cancel interference from sidelobes. For instance, creating a secondary array 
to preferentially receive an offending signal from one of the main array’s 
sidelobes with the intent of subtracting this signal from the main array’s 
output. Other applications of multiple arrays include more advanced 
direction of arrival estimation, including triangulating separate array’s 
estimates to yield a direction and distance (which can be tracked over time 
to get a direction of motion estimate). 


Code and Poster 
A link to our MATLAB code and a PDF version of our poster 


Code and Poster 


For those who are curious, we’ve provided a link to our MATLAB code and 
poster. Both of them can be found here. 


